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Abstract 

Defect-free and defected composite thin shells with ply orientation (90/0/ ± 75) made of graphite /epoxy are simulated 
for damage progression and fracture due to internal pressure and axial loading. The thin shells have a cylindrical geometry 
with one end fixed and the other free. The applied load consists of an internal pressure in conjunction with an axial load at 
the free end, the cure temperature was 177°C (350°F) and the operational temperature was 21°C (70°F). The residual stresses 
due to the processing are taken into account. Shells with defect and without defects were examined by using CODSTRAN an 
integrated computer code that couples composite mechanics, finite element and account for all possible failure modes 
inherent in composites. CODSTRAN traces damage initiation, growth, accumulation, damage propagation and the final 
fracture of the structure. The results show that damage initiation started with matrix failure while damage /fracture 
progression occurred due to additional matrix failure and fiber fracture. The burst pressure of the (90/0/ + 75) defected 
shell was 0.092% of that of the free defect. Finally the results of the damage progression of the (90/0/ ± 75), defective 
composite shell was compared with the (90/0/± 0), where 0= 45 and 60, layup configurations. It was shown that the 
examined laminate (90/0/ ± 75) has the least damage tolerant of the two compared defective shells with the (90/0/ + 0), 
0= 45 and 60 laminates. 


1. Introduction 

Aircraft, marine and automotive vehicle industries 
use composite shells because of their low weight and 
high stiffness and stability features. Design consider- 
ations with regard to the durability of composite 
shells require a priori evaluation of damage initiation 
and propagation mechanisms under expected service 
loads. Concerns for safety and survivability of criti- 
cal components require quantification of the compos- 
ite structural damage tolerance during overloads. 


Corresponding author. Fax: 4- 1-216-4333252. 


Characteristic flexibilities in the tailoring of compos- 
ite structures make them more versatile for fulfilling 
structural design requirements. However, these same 
design flexibilities render the assessment of compos- 
ite structural response and durability more complex, 
prolonging the design and certification process and 
adding to the cost of the final product. It is difficult 
to evaluate composite structures because of the com- 
plexities in predicting their overall congruity and 
performance, especially when structural degradation 
and damage propagation take place. The predictions 
of damage initiation, damage growth, and propaga- 
tion to fracture are important in evaluating the load 
carrying capacity, damage tolerance, safety, and reli- 
ability of composite structures. The most effective 
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way to obtain this quantification is through inte- 
grated computer codes that couple composite me- 
chanics with structural analysis and damage progres- 
sion models. The Composite Durability STRuctural 
ANalysis (CODSTRAN) computer code has been 
developed for this purpose by integrating and cou- 
pling the following disciplines: (i) mechanics of 
composites, (ii) structural analysis (FEM) and (iii) 
damage progression tracking. CODSTRAN computer 
code was used to simulate the damage progression in 
a variety of fiber composite structures such as: pro- 
gressive fracture of fiber composite thin shell struc- 
tures [1-3], fiber composite stiffshed panels [4], 
Stiflened adhesively bonded composite structures [5], 
dynamic damage progression of a containment struc- 
ture hit by an escaping blade [6], damage progression 
in adhesively bonded fiber composite thin shell 
structures [7] and [8], damage progression in bolted 
composite structures [9] and simulation of the 
Iosipescu shear testing [10]. 

The objective of the present paper is to demon- 
strate what can be accomplished by integrating com- 
posite mechanics, finite element structural analysis, 
and tracking of composite failure modes into a 
stand-alone computer code. The computer code in 
our case is CODSTRAN. Reference to CODSTRAN 
throughout this report should not be construed nei- 
ther as NASA endorsement nor as an advertisement 
of that computer code. 


2. Progressive fracture methodology in codstran 

CODSTRAN is an integrated computer code 
which was developed by coupling three modules: 
composite mechanics (ICAN [11]), finite element 
analysis (MHOST [12]) and a damage progression 
modelling algorithm: 

(1) ICAN is a composite mechanics computer 
code [11] that provides the constituent (fiber and 
matrix) material properties using an available data 
bank, and computes the ply properties and the com- 
posite properties (effective properties) of the lami- 
nate in a hygrothermal environment. The theory of 
the code is based on the micromechanics of compos- 
ites and the laminate theory. ICAN has the ability to 
compute the ply stresses by knowing the stress resul- 
tants (force per laminate thickness, where force can 


PLY FAILURE CRITERIA 
♦ MAXIMUM STRESS FAILURE CRITERION. 


THE SIX PLY STRESS COMPONENTS ARE ALONG THE MATERIAL AXES . 

Sluc < °lh * Slut 

S L22C < °L22 < S L22T 
S L33C < °Li3 < S L33T 
S L12<-) < °LU « ®L12<+) 

S L23{-) < a L23 * Sl23<+) 

S LU(-) < °LU < S L13<*) 



♦ MDE COMBINED STRESS FAILURE CRITERION. 

F*l-[ (olhj/Slh*) 2 + (o L22b /S L22l>> 2 ~ K LlW a Lllii /S Lllj . ] 

+ ^°LI2S^L12S^ *• 

Fig. 1. Ply failure criteria. 


be a concentrated load, a bending or a twisting load). 
In ICAN failure criteria were established (Fig. 1), for 
the detection of the ply failures as follows: (a) the 
maximum stress criterion, in which failure occurs 
when the individual ply stress a Lij for ij= 1,2,3, 
exceeds the respective ply strength S LiJ for ij = 
1,2,3; and (b) the modified distortion energy crite- 
rion, in which the combination of the ply stresses is 
taken into account. In Fig. 1, a and b are referred to 
the tensile and compressive stresses, respectively. In 
both criteria the ply stresses are referred to the 
material axes 1, 2, 3, and the direction of the 0° 
fibers are along the direction of the material 1-axis. 
For example a laminate with configuration 
(90/0/ ±75) and ply stresses at the top ply (90°) 
are shown in Fig. 2. In ICAN, the described failure 
modes of the plies are: failure due to the fiber 


Pi; strewcs at the top ply (90°) 



°L22 * Ply Iruuvaic stren 
*L12 • Ply their Ureu 

Fig. 2. A typical laminate configuration (90/0/ ±75) and the ply 
stresses at the top ply (90°). 
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Fig. 3. CODSTRAN progressive fracture simulation cycle. 


fracture in tension or in compression; damage due to 
the matrix fracture in tension or in compression; and 
damage due to intralaminar and interlaminar shear 
fracture. 

(2) MHOST is a finite element computer code 
[12] for the solution of structural analysis problems. 
The code has the capability to perform linear or 
nonlinear static and dynamic analysis. MHOST has a 
library with a variety of elements and for the present 
work the four node shell element was used. By 
supplying the boundary conditions, the desire type of 
analysis, the applied loads and the laminate proper- 
ties (using ICAN) MHOST performs the structural 
analysis. In addition MHOST provides the computed 
stress resultants to the ICAN code; and then ICAN 
computes the developed ply stresses for each ply and 
checks for ply failure. 

(3) A module that keeps track the failure modes 
and communicate these modes to ICAN in order to 
degrade the properties associated with the respected 
failure modes. 

An integrated schematic of the CODSTRAN sim- 
ulation cycle is shown in Fig. 3. In this figure, from 
the left side along the clockwise direction the mate- 
rial properties of the constituents (Fiber and matrix) 
are provided by ICAN’s data bank, next the ply 
properties are computed by using the micromechan- 


ics theory, and the laminate properties are computed 
using the laminate theory. These properties in con- 
junction with the finite element mesh, the loads and 
the boundary conditions are incorporated into 
MHOST. MHOST performs the structural analysis 
and provides the computed stress resultants in ICAN 
(in the right side of Fig. 3), where ICAN proceeds to 
compute the ply stresses using the laminate theory 
and checks for ply failure. 

The nonlinear structural analysis in MHOST code 
is performed in conjunction with an incremental load 
algorithm. The load is increased in small increments 


CODSTRAN LOAD INCREMENTATION 



Fig. 4. CODSTRAN load increment. 
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(equilibrium positions). Within each equilibrium po- 
sition a number of iterations are performed for dam- 
age growth and accompanying internal structural re- 
sponse to sustain the applied load. Fig. 4. In each 
iteration, the structure is checked for ply failure. If 
damages are detected in the structure, the model is 
automatically updated with a new finite element 
mesh and new laminate properties, and a new finite 
element analysis is performed; the above iterations 
continue until no further damage occurs (equilibrium 
position). After that, the load is increased, and the 
above procedure is repeated until the final failure of 
the structure. Following the above procedure, the 
damage progression, fracture and collapse of the 
structure is monitored. 


3. Geometry, loads, materials and the FE model 
of the thin shell 

The geometry of the fiber composite shell is a 
thin cylinder. The thin cylinder has inner diameter of 
305 mm (12 in.) and length of 760 mm (30 in.). The 
boundary conditions are one end fixed supported and 
the other free. The corresponding graphite fibers 
(AS-4) and epoxy matrix (HMES) properties of the 
fiber composite thin shell were obtained from a data 

Table 1 

AS-4 graphite fiber properties 

Number of fibers per end = 10000 
Fiber diameter = 0.00762 mm (0.3E-3 in) 

Fiber density = 4.04 X E — 7 kg/m 3 (0.063 lb/in 3 ) 

Longitudinal normal modulus — 227 GPa (32.9 X E 4- 6 psi) 

Transverse normal modulus = 1 3.7 GPa ( 1 .99 X E + 6 psi) 

Poisson’s ratio (y, 2 ) = 0.2 

Poisson’s ratio (r^j) = 0.25 

Shear modulus (G 12 ) = 13.8 GPa (2XE + 6 psi) 

Shear modulus (C 2? ) = 6.9 GPa (1 XE + 6 psi) 

Longitudinal thermal expansion coefficient = 1 XE-6/°C 
( - 0.55 X E - 6/°F) 

Transverse thermal expansion coefficient = 1 XE — 6/°C 
(-0.56XE-6/°F) 

Longitudinal heat conductivity = 43.4J m/h/m 2 /°C 
(580 BTU-in/h/in 2 /°F) 

Transverse heat conductivity = 4.34 J m/h/m 2 /°C 
(58 BTU-in/h/in 2 /°F) 

Heat capacity = 712 J/kg/°C (0.17 BTU/lb/°F) 

Tensile strength = 3723 MPa (540 psi) 

Compressive strength = 3351 MPa (486 psi) 


Table 2 

HMHS epoxy matrix properties 

Matrix density = 3.4XE-7 kg/m 3 (0.0457 lb/in. 3 ) 

Normal modulus = 4.27 GPa (629 psi) 

Poisson’s ratio = 0.34 

Coefficient of thermal expansion = 0.72/°C (0.4 X E - 4/°F) 
Heat conductivity = 1.25 BTU-in./h/in. 2 /°F 
Heat capacity = 0.25 BTU/lb/°F 
Tensile strength = 84.8 MPa (12.3 psi) 

Compressive strength — 423 MPa (61.3 psi) 

Shear strength = 1 48 MPa (2 1 .4 psi) 

Allowable tensile strain = 0.02 
Allowable compressive strain = 0.05 
Allowable shear strain = 0.04 
Allowable torsional strain = 0.04 
Void conductivity = 16.8 J m/h/m 2 /°C 
(0.225 BTU-in./h /in. 2 /°F) 

Glass transition temperature - 216°C (420 °F) 


bank of composite constituent material properties 
resident in CODSTRAN and are given in Tables 1 
and 2, respectively. The AS4/HMES ply strengths 
computed by ICAN code are given in Table 3. The 
fiber composite thin shell consists of eight 0.136 mm 
(0.00535 in.) plies resulting in a composite shell 
thickness of 1,088 mm (0.0428 in.). The laminate 
configuration is (90/0/ + 75), with the 90° plies in 
the hoop (or circumferential) direction and the 0° 
plies in the axial direction. Fiber orientations in the 
±75° is shown in Fig. 2. The fiber volume ratio is 
60% and the void volume ratio is 2% of the total 
volume of the structure. The residual stresses is been 
taken into account due to the processing. The cure 
temperature was 177°C (350°F) and the pressuriza- 
tion/use temperature is 21°C (70°F). The moisture 
content was zero. The closed-end cylindrical pres- 
sure vessel is simulated by applying a uniformly 

Table 3 

AS-4/RMHS ply strengths 

S u 1T = 1930.30 MPa (280 psi) 

S L j ic ~ 1475.85 MPa (210 psi) 

S L22T — 91.38 Mpa (13 psi) 

5 L22C = 228 27 MPa ( 33 P si ) 

S u 2 -65.57 MPa (9.5 psi) 

^3 = 59.98 MPa (8.7 psi) 

Where 1 , 2, 3 are the material axes of the ply. The direction of the 
fibers are parallel to 1-axis. T is for tension and C is for 
compression. 
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FIBER COMPOSITE THIN CYLINDRICAL STRUCTURE 

CYLINDER LENGTH, L = 76 J cm (30 in.) 

INNER RADIUS, r = 15.24 cm (6 in.) 

LAMINATE THICKNESS, t = 0.11 cm (0.043 in.) 

GRAPHITE/ EPOXY , (90/0/175)^ FIBER VOLUME FRACTION = 60 % 



Fig. 5. Applied loads and boundary conditions of the thin shell structures. 


distributed axial tension such that the generalized 
axial stresses in the shell wall were half of those 
developed in the hoop direction, Fig. 5. The shell is 
subjected to a monotonically increasing internal 
pressure until it burst. Shells with defect and defect 
free are examined. The defective shell has a through- 
the-thickness defect at a node located at the half 
length of the shell. After the load is applied in the 


defective shell, a crack is formed at the location of 
the defective node as shown in Fig. 5. The geometry 
of the formed crack consists of a 95 mm (3.75 in.) 
long thin axial slit that is superimposed on a 60 mm 
(2.36 in.) long circumferential slit. Each finite ele- 
ment model contained 544 nodes and 512 uniformly 
sized rectangular elements. Fig. 6. Damage initiation, 
growth, fracture progression, and global structural 


DAMAGE INITIATION , AT INTERNAL PRESSURE 25 PSI. 



GRAPHITE /EPOXY. (90/0/75/-75)s 


Fig. 6. Finite element mesh. 
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fracture stages were investigated. Computed results 
are presented up to global fracture for defect-free 
shells and shells with through-the-thickness defects. 


4. Results and discussion 

(1) Fiber composite (90/0/ ±75), free defects 
thin shells. The damage of the thin shell structure is 
plotted versus the internal pressure in Fig. 7. 1 . Dam- 
age is defined as the ratio of the volume of the 
damaged plies divided by the total volume of the 


plies. The depicted points of Fig. 7.1 are equilibrium 
points that were discussed in the previous chapter in 
CODSTRAN Methodology. Damage initiation oc- 
curred at pressure of 1.06 MPa (153 psi) with matrix 
cracking at the outermost ply 2 (0°). The damage 
progressed with shear failure at the ply 3 (75°) and 
matrix cracking at the plies 7 (0°) and 8 (90°). 
Further, the damage progressed and matrix cracking 
occurred at the outermost ply 1 (90°), and at the 5 
( — 75°) and 6 (75°). After the initial matrix degrada- 
tion stage the pressure was increased without activat- 
ing additional damage modes. The thin shell burst 


(a) 


DEFECT FREE THIN CYLINDER (90/0/±75) t 


10 r 


0.8 


0.6 

uj 

O 

< 

5 

< 

Q 04 


0,2 


0.0 



DAMAGE INITIATION 

PLY 2(0°) MATRIX FAILURE^ L22T(lst) 
PLY 3 (75®) SHEAR FAILURE,a U2S (2n 
PLY 7(0°) MATRIX FAILURE,o U2T 
PLY 8 (90°) MATRIX FAILURE^.™ 
SHEAR FAILURE^ L ^ 

PLY 1(90°) MATRIX FAILURE, o U2T(3rd) 
PLY 5 (-75°) MATRIX FAILURE^ L22T 
PLY 6 (7 5°) MATRIX FAlLURE,a L22T 


3 4 5 

INTERNAL PRESSURE (MPa) 


(b) 



Normalized pressure (Max. pressures .367 ksi) 


Fig. 7. (1) Damage versus pressure in defect-free thin shell. (2) effect of the laminate configuration, on the damage of the defect-free thin 
shells. 
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Fig. 8. ( 1 ) Damage versus pressure in defective thin shell. (2) Effect of the laminate configuration, on the damage of the defected thin shells. 


when ply 2, the outermost 0° axial ply experienced 
fiber fractures that caused rapid damage propagation 
to the other plies and structural fracture. Hie struc- 
tural fracture was initiated by the tensile failure of 
axial fibers. The simulated burst pressure was 7.83 
MPa (1.14 psi). 

The damage progression of the (90/0/ ±75) 
laminate compared with the (90/0/ ±45), and 
(90/0/ ±60), layup configurations in Fig. 7.2. In 
the case of the (90/0/ ±45), laminate, damage 
initiation occurred at 0.8 MPa (116 psi) with matrix 
cracking in the first ply (0°) followed by the 45° ply 
immediately adjacent to it, then the innermost (0°) 


ply and the 90° ply on the interior surface. After 
damage initiation, all plies gradually sustained ma- 
trix cracking as the pressure was increased to 0.84 
MPa (121 psi). After this damage accumulation stage 
the pressure was significantly increased without acti- 
vating additional damage modes. The thin shell burst 
when the 90 hoop plies experience fiber fractures at 
4.8 MPa (696 psi). In the case of the (90/0/ ± 60), 
laminate, initial damage occurred at 0.97 MPa (140 
psi) with matrix cracking at the 2nd ply (0) followed 
by the 3rd ply (60) then by the 7th ply (0) and finally 
by the 8th ply (90). After damage initiation, all plies 
gradually sustained matrix cracking as the pressure 


DAMAGE INITIATION , AT INTERNAL PRESSURE 25 PSI. 




GRAPHITE /EPOXY (90/0/75/-75)s 


Fig. 9. Damage initiation at 0.17 MPa (25 psi). 
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DAMAGE PROPAGATION , AT INTERNAL PRESSURE 37.5 PSI 




GRAPHITE/EPOXY , (90/0/75/-75)s 


Fig. 10. Damage progression at 0.26 MPa (37.5 psi). 


was retained virtually constant. After this damage 
accumulation stage, the pressure was significantly 
increased without activating additional damage 
modes. The thin shell burst when the outermost 90° 
hoop plies experienced fiber fractures at 9.42 MPa 
(1366 psi). 

(2) (90/0/ ±75), fiber composite defective thin 
shells. The damage versus the internal pressure is 
shown in Fig. 8.1. The damage initiated at internal 
pressure of 0.172 MPa (25 psi.) with matrix fracture 


along the hoop direction nodes 272 and 274 in the 
7th (0°) ply and along the axial direction nodes 241 
and 305 in the 6th (75°) and 8th (90°) plies. The 
damage progressed continuously with shear failure at 
the 3rd (75°) ply in both hoop and axial nodes, and 
with matrix failure at the 8th (90°) ply at the hoop 
and at the 5th (-75°) ply, at the axial node. The 
damage of the plies continued with relatively small 
increases in pressure in the hoop direction until fiber 
fracture occurred at all plies except of the 5th ( — 75°). 


DAMAGE PROPAGATION, AT INTERNAL PRESSURE 62.5 PSI 



i 

iz 


Fig. 1 1. Damage progression at 0.43 MPa (62.5 psi). 
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DAMAGE PROPAGATION, AT INTERNAL PRESSURE 70.8 PSI 



!y 

\ 

iz 


~.x 


GRAPHITE /EPOXY 


(90/0/75/-75)s 


Fig. !2. Damage progression at 0.49 MPa (70.8 psi). 


Increasing the load further, the damage of the shell 
progressed steadily until it reached the value of 0.62 
MPa (90 psi) where the fracture growth became 
unstable, increased rapidly, and structural fracture 
occurred at burst pressure equal to 0.717 MPa (104 
psi). 

The damage progression of the (90/0/ ±75) 
laminate compared with the (90/0/ ±45) and 
(90/0/ ± 60) layup configurations as shown in Fig. 
8.2. In the case of (90/0/ ± 45) laminate, the dam- 


age initiation started at pressure 0.175 MPa (25 psi). 
The damage growth localized adjacent to die defect 
until a 0.74 MPa (108 psi) pressure reached. Damage 
propagation by fiber fractures concentrated near the 
defect. The thin shell burst at pressure equal to 1.58 
MPa (230 psi). In the case of (90/0/ ±60) lami- 
nate, damage initiation started at pressure 0.17 MPa 
(25 psi) by matrix cracking adjacent to the defect at 
the 8th ply (90°) at the circumferential slit tips and 
by the 7th ply (0°) at the axial slit tips. Damage 


DAMAGE PROPAGATION , AT INTERNAL PRESSURE 90. PSI 



Q RAP HITE/E POXY, (90/0/75/-75)s 


k 


Fig. 1 3. Damage progression at 0.62 MPa (90.0 psi). 
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growth remained localized adjacent to the defect as 
the pressure was increased to 0.91 MPa (132 psi). 
The thin shell burst at 1.9 MPa (275 psi). 

Regarding the comparison of the direction of the 
damage progression: the damage growth of the 
(90/0/ ± 45) laminate was mainly in the axial direc- 


tion of the shell. The damage growth of the 
(90/0/ ± 75) laminate was mainly in the circumfer- 
ential direction of the shell. Finally, the direction for 
the damage growth for the (90/0/ ± 60) laminate 
was both in the circumferential as well as in the axial 
direction. 


TOP PLY (90) PLY LONGITUDINAL STRESSES 
DAMAGE INITIATION , AT INTERNAL PRESSURE 25 PSI. 



GRAPHITE /EPOXY. 



(90/0/75/-75)s 
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Fig. 15. At the top ply (90°), ply longitudinal stresses at 0.17 MPa (25 psi). 
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TOP PLY (90) PLY LONGITUDINAL STRESSES 
DAMAGE PROGRESSION AT INTERNAL PRESSURE 70.8 PSI 
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Fig. 16. At the top ply (90°), ply longitudinal stresses at 0.49 MPa (70.8 psi). 
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Fig. 17. At the top ply (90°). ply longitudinal stresses at 0.62 MPa (90 psi). 
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TOP PLY (90) PLY TRANSVERSE STRESSES 
DAMAGE PROGRESSION AT INTERNAL PRESSURE 90. PSI 
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Fig. 20. At the top ply (90°), ply transverse stresses at 0.62 MPa (90 psi). 


PSI 

15845 . 

1 4247 . 
12650 . 
11053 . 
9456 . 
7859 . 
6262 
4665 
3068 . 
1470 . 
- 126.8 
- 1724 . 
-3321 
- 4918 . 
- 6515 . 
- 8112 . 


The fracture progression of the (90/0/ ± 75) thin 
shell structure is depicted in Figs. 9-14, for different 
load increments. At the top ply (90°), the ply longitu- 
dinal stresses (in the hoop direction) are plotted in 
Figs. 15-17; and the ply transverse stresses (in the 
axial direction) are plotted in Figs. 18-20 for differ- 
ent load increments. At the damage initiation stage 
high ply stresses were developed at the edges of the 
crack, and when the fracture progressed the distribu- 
tion of the ply stresses around the fractured area 
changed and decreased. 


5. Summary 

In this investigation, the simulation of the struc- 
tural and damage progression response of a 
(90/0/ ±75), composite thin shell structure was 
examined. The damage progression of the (90/0/ + 
75) laminate was compared with the (90/0/ ± 6 ), 
(for 0 — 45 and 60) layup configurations. The signif- 
icant results are as follows: 


(1) For the case of the defect free (90/0/ ± 75) 
thin shell structure damage initiation started at inter- 
nal pressure 0.14% of the burst pressure 7.83 MPa 
(1.1 psi) with matrix cracking at the outermost ply 2 
(0°) and followed by the same type of failure at the 
plies 3 (75°), 7 (0°) and 8 (90°). The damage pro- 
gressed and when fiber fractured occurred at the 
outermost ply 2 (0°) rapid fracture initiated and 
collapse of the structure occurred at 0.717 MPa 
(0.104 psi). 

(2) For the case of the defected (90/0/ ± 75) thin 
shell structure, the damage initiation started at inter- 
nal pressure of 0.24% of the burst pressure 0.717 
MPa (0.104 psi), with matrix cracking at the plies 
near the nodes of the crack tips (hoop and the axial 
directions nodes). Increasing the applied pressure, 
the fiber fracture progressed in the hoop direction at 
a faster rate than in the axial direction, until the 
collapse of the structure. 

(3) Comparing the defected with the defect free 
thin shell structure, it was observed that damage 
initiation started at 0.16% of the defect free thin shell 
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pressure 1.06 MPa (0.153 psi), while burst pressure 
occurred at 0.991% of the defect free thin shell 
pressure 7.83 MPa (1.13 psi). 

(4) The direction of damage growth for the 
(90/0/ + 75) laminate was mainly in the circumfer- 
ential direction of the defected shell. This laminate 
was the least damage tolerant of the two compared 
defective shells with (90/0/ ± 45) and (90/0/ ± 60) 
layup configurations. 

(5) The direction of damage growth for the 
(90/0/ ± 45) laminate was mainly in the axial direc- 
tion of the defective shell, had significant damage 
tolerance and sustained the maximum amount of 
damage prior to bursting. 

(6) The direction of damage growth for the 
(90/0/ ± 60) laminate was both in the circumferen- 
tial as well as in the axial direction of the defected 
shell. This laminate was able to sustain the highest 
burst pressure. However, the amount of damage 
sustained was lower than that corresponding to the 
(90/0/ + 45) laminate. 
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